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Abstract: The quantum evolution of a model of the universe with account of two scalar 
fields (dilaton and inflaton) is considered. For this case, the closed and flat models has 
been examined. It is shown that in both cases the realization of conditions necessary for 
inflation is strongly depends on coupling constant for dilaton and inflaton fields. 
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1. Introduction 

In Ref. (l| the new inflationary scenario (so-called extended inflation) based on a combina- 
tion of old inflation Q with Jordan-Brans-Dicke gravitational theory || (hereafter JBD) 
has been offered. It allows to slow down a speed of expansion of the universe at the initial 
stage of its evolution. Thus there was power-law expansion of the universe (instead of 
exponential one). Such delay has allowed to solve a number of problems of the old inflation 
noted still in Refs. |fflj and Q (in particular, the problem of graceful exit). Unfortunately 
in model of La and Steinhardt there was the basic difficulty connected with impossibility of 
an establishment of thermal equilibrium in walls of bubbles for cosmologically reasonable 
times. It results in violation of isotropy of background radiation. With the purpose of 
avoidance of this problem, it is necessary that the coupling constant in JBD theory has 
been limited from above uj < 20 |J. However, the lower observable bound is uj > 500 §. 
For resolving of this problem, Damour, Gibbons and Gundlach Q (hereafter DGG) has 
proposed the model in which, except of visible matter, the invisible (field) component is 
present as well. Action in their model has a form (in Jordan (physical) frame): 

SIsV'^'Vv,^/] = S B D[g f w,$] +S v \g^,ipv] + Si[g IJ , u ,^,ifji], (1.1) 
where <& is JBD scalar field (dilaton), Sbd[9/iu, 3>] - Brans-Dicke action || 



1.2) 



and Sy, Sj are actions for visible tpy and invisible ipj components of material fields. As 
the usual matter is described as an ideal fluid, it does not play any role during inflation at 
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this stage of evolution of the universe and it is possible to dropped it. Invisible component 
in DGG model reads: 

Silg^Q^i] = J d 4 x^ [^-V^Vj^j - ^-flVtyj)] > (1-3) 

where (3 = fii/fiy and (3yj are coupling constants of visible and invisible matter with 
dilaton field. (Hereafter we work in geometrical units, i.e. c = 1 and 8ttGn = 1, where 
Gn is the current value of Newtonian constant.) DGG worked in so-called Einstein frame 
which turns out from Eq. ( |1.1[ ) via conformal transformation of the metrics g^ u = 2Qg^ u 
with simultaneous redefinition of the variable $ = | exp ((f)/ fa), where 0o = \/^> + 3/2. 
Using Eq. (1.1) Holman, Kolb and Wang B has shown that at a choice of inflaton field as 



an invisible matter, it is possible to find region on a plane u — (3 in which the mentioned 
above restrictions uj < 20 (necessary for successful inflation) and u> > 500 (claimed from 
observations) can be satisfied simultaneously. It is possible because, instead of standard 
restriction u < 20 from model Q, the new relation to/ f3 2 < 20 will work. It can be easily 
satisfied at u > 500. 

The further development of the model was reduced to introduction in Eq. (|1.3j ) arbi- 
trary exponents n and m instead of (1 — (3) and 2(1-/3) || that has allowed to investigate 
the generalized process of nucleation and percolation and to calculate time-depended nu- 
cleation rate per unit volume. In Ref. |l(]] the restrictions on parameters n and m at which 
the sufficient inflation was provided have been found. 

It seems interesting to consider quantum evolution for the model of extended inflation. 
The main purpose of this paper is to find necessary conditions for ensuring of inflation on 
post-quantum stage. In Section 2 the Wheeler-DeWitt equation for closed, open and flat 
model is obtained in general form. Section 3 is devoted to problem statement and to general 
description of solution technique. The derivation and solving of corresponding equations 
for closed and flat models are presented in Section 4. Finally, in Section 5 the probabilities 
of tunneling and above the barrier reflection for closed and flat models accordingly are 
found. 

2. Wheeler-DeWitt Equation 

We will work in Einstein frame. Thus, using conformal transformation mentioned in In- 
troduction and using designation b = 1 — @, we can obtain from Eqs. ( |1.2[ ) and ( |Q| ) the 
following Lagrangian (see also Ref. Ji~l|]): 



£ = \>fl[-R + gTd^t + ^gTd^a - e 2b ^m 2 a 2 }. (2.1) 

Here the massive inflaton field a with the simplest form of potential V = \m 2 a 2 is using 
as invisible matter. Let us write out the general expression for metrics for closed, flat and 
open universes: 

2 



ds 2 = N 2 (t)dt 2 - a 2 (t) 



(2.2) 
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In this expression k = —1,0,1 are for closed, flat and open models accordingly, N(t) is 
lapse function. As is known, the metrics in itself does not yet defines global topology of 
space-time. It enables to consider spacelike sections as compact in a case of open and flat 
models. Thus the volume, limited by spacelike hypersurfaces, remains finite and have a 
form: 



a! 3 an 



-h. 



(2.3) 



(h is the three-dimensional metrics.) In a case of homogeneous model, using Eq. fl2,2|) , we 
can obtain: 



V k = a v k , v k 



sin* 



k 



■ sm ( 



(2.4) 



where the value of v k depends on topology of the spacelike sections (see, e.g., Ref. p2| 



and references inside). Then, using Eqs. (2.1) and ( |2.2| ), we have the following Lagrange 
function: 



L = 3v k N 



1 2 

ka - j^aa 



1 



1 



ob<f> 



Further, using Eq. ( |2.5[ ) and calculating the conjugate momenta 

Qv k L v k 3 



Pa 



■— aa, P^ = j^a 



Per = ^a 3 e b % 



we can obtain the corresponding canonical Hamiltonian: 

1 

2"^ 3 



A = ™ 

Vk 



3vla (k- im 2 a W 
6 



(2.5) 



(2.6) 



(2.7) 



Quantization of Eq. Q2.7|) is reduced to replacement of the conjugate momenta p a ,P<f>,Pa 
on operators —id /da, —id/d(j) and —id /da accordingly. Wheeler-DeWitt equation turns 
out at action of the Hamilton operator on the wave function ifj(a, (f>, a) and looks like: 



12 da \ da + 2a 2 



+ e~ 



^ + 3v 2 k a 2 (k - l -m 2 a 2 e 2b ^^j = 0, (2. 



d 



where p is factor ordering which introduction is caused by ambiguity of commutation 
properties of the operators a and p a . Transformation of wave function 

allows us to exclude the first derivative in Eq. ([^) that gives: 



1 d 2 ip 1 (d 2 i\) 



+ 



12 da 2 la 2 



+ e~ 



,d 2 ^ 
da 2 



+ 



24a 2 



P 



+ 3via 2 ( k 



6 



:m 2 aW 



ip = 0. 
(2.9) 

Let us choose the factor ordering p = 1. Also we shall use the following redefinitions: 
(/>, cr — > \/6<p,a, b — > y/l/6b, v k — > l/6vk- Finally, from Eq. fl2.9|) we have: 



d 2 V> 1 /^5 2 ^ 
9a 2 a 2 



-6(6 



d 2 v> 
d^ 2 



+ = 0, 



where superpotential 



IF 



1 

4a 2 



2 2 / i 2 26(4 2 2 

+ v£a ( k — a e v m a 



(2.10) 
(2.11) 



- 3 - 



3. Problem statement 



The obtained equation describes quantum evolution of the universe. Thus the statement of 
problem essentially depends on type of a model of the universe which is examined: closed, 
flat or open. In view of equality to zero of full energy of the universe, it is clear from form 
of superpotential ( 2.11 ) that consideration of the problem is reduced to research either the 
process of tunneling in case of the closed model or the process of above the barrier reflection 
of the universe's wave function in case of flat and open models. Such consideration was 
carried out repeatedly in various statements. However, in most cases one-dimensional 
problems were examined. In this cases wave function depends only on one variable - scale 
factor. At that it was possible to obtain exact solutions in a number of problems |13| . 
And, that is even more important, in one-dimensional statement it is seems possible to 
find clear physical sense of the wave function. In a case of multidimensional problems, 
when the wave function depends also from field variables, the problem is complicated not 
only by mathematical difficulties but also by complexities with understanding of physical 
sense of the multidimensional wave function. In this connection some authors proposed 
the methods of reducing of multidimensional problems to one-dimensional one. It allows 
to give to the wave function physically clear sense fl4]| . 

In the given work we shall use one of such ways which was suggested in Ref. [15]. 
The essence of this method consists in consideration of semi-classical evolution of a model 
of the universe. In this case the corresponding Hamilton-Jacobi equation, similar to the 
equation from classical mechanics, is written down for the Wheeler-DeWitt equation. This 
equation is nonlinear first-order partial differential equation. For its solution the known 
method is used. The essence of this method consists in representation of the given non- 
linear partial differential equation through a system of the first-order ordinary differential 
equations (so-called set of characteristic equations). The given equations describe a set 
of characteristics which are form an integral surface of initial partial differential equation. 
Distinctive feature of the system of characteristic equations is that the functions, entering 
in equations, depend only on one argument. Thus, initially multidimensional problem is 
reduced to one-dimensional one that allows to find physically sensible one-dimensional wave 
functions along each of characteristics. Then for calculation, for example, of permeability 
of a barrier, it is necessary only to sum and to average permeability on all characteristics. 



4. Derivation and solving of a set of characteristics 

Let us pass directly to research of our model. 
Closed model 

At the beginning, let us consider the case of closed universe. As it was already specified 
above, in this case the process of tunneling of the wave function will be take place. In 
semi-classical treatment of the given problem we can present the wave function as t/j = e~ s . 
Then, being limited to zero approximation, we obtain the following equation for action S 
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from Eq. (|2.10[) : 



dS_ 

da 



or 



dsy 



M (dS 
\da 



+ W = 0. 



(4.1) 



The corresponding set of characteristics under a barrier will be (see., e.g., Ref. fllql ): 



d<j) 
da 
dq 

da 
dS_ 

da 



a 2 F 

1 be 

~2 
W 

T 



-b,p 



da 
da 

,2 



a?F 



a 2 F 



bv 2 a 4 e 2b *m 2 a 2 
F 



dv 
da 



v 2 a 4 e 2b ^m 2 a 
F 



(4.2) 



where q = dS/ dcj) and v = dS/ da are the generalized momenta and F = \J (q 2 + v 2 )/a 2 + W. 

This equations describe a process of tunneling of a "particle" along the characteris- 
tics. Each of such characteristics differs from another by definition of corresponding initial 
conditions. For this purpose let us analyse the potential ( |2.11| ) at k = 1. As we see, the 
whole minisuperspace is divided into three areas (see Fig. @): 1) internal area, where 



W(a)k 




Figure 1: The conventional form of an one-dimensional barrier is shown. Tunneling of the wave 
function occurs from the point ao, where initial conditions are set, into the point of an output from 
under the barrier a±. The initial conditions in ao define the form of a barrier which is individual 
for each characteristic. 

there is no classical space-time and the metrics is under the strong quantum fluctuations; 
2) classically forbidden area under a barrier; 3) classically allowed area. Initial conditions 
are set on bound between internal and classically forbidden areas. For finding of the bound, 



it is necessary to equate superpotential ( 2.11 ) to zero. In generally, it will be a surface. 
However, if to choose as initial conditions (p,a = 0, evolution of the model will begin in a 
point ao = (l/4v^) 1 / 4 , where v\ is a three-volume of the closed universe. In case of sphere 
v\ = 2tt 2 . Further, we should set initial momenta q = qo and v = vq which will vary from 
up to oo. So, all characteristics will be differ by different values of momenta. Using 
the specified initial conditions and solving numerically the set of equations ( [4.2|) we have 
obtained the following results (see Figs. (||)). 
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Figure 2: The dependence 4>(a) (left panel) and a(a) (right panel) for different b: 1 - for b = 1. 
2 - for b = 0.2, 3 - for b = 0.02. Dotted lines are pointing out the points of output from under a 
barrier. 



This implyies that in process of under a barrier evolution, there is a growth both 
inflaton and dilaton fields. Thus it is important to note that thickness of a barrier and, 
consequently, the rate of growth of fields strongly depends on parameter b: the less value 
of b means the more width of a barrier and the more value of fields on an output from 
under a barrier. On the other hand it is clear that after an output from under a barrier the 
parameters of the model should provide conditions for the subsequent inflationary stage. 
At least, there are two such conditions: 1) inflaton field should be the Planck's order; 
2) it should varies slowly enough to ensure long inflationary stage. In the given model 
such conditions are provided in a case when the coupling constant b is small enough. This 
conclusion is in agreement with specified in Ref. pi restriction on value of b << 1. 



Flat model 

It is possible also to consider evolution of the given model for cases when k = 0, — 1. It is 
obvious that in these cases the potential barrier will be lay in negative area and only above 
the barrier reflection is possible. In view of absence of principle distinctions of process of 
above the barrier reflection for open and flat model, we shall consider, for example, an 
evolution of flat model (i.e. we should put in Eq. (ETUI) k = 0). Since full energy of the 
universe is equal to zero and the superpotential always is less than zero, then in case of 
flat model evolution of the universe will always take place in classically allowed area. In 
this connection let us search for solution of equation (|2.9|) as ip = e~ lS . Then we obtain: 



da ) 



1 +h 



dS\ 



+ e 



W = 0. 



(4.3) 



Then the set of characteristics reads: 
deft q 



da 



a 2 F : 



da 
da 



a 2 F 
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dq 1 be- b<t> v 2 ba 4 e 2b ^m 2 a 2 dv _ a 4 e 2b( ^m 2 a 

~da~~2 a 2 F + F ' F ' 

dS_ _ _W_ 

da~ ~ ~F' 

where F = sj (q 2 + v 2 )/a 2 — W . (Here we have choose Vq = 1 (see Ref. Jl2[).) For 
numerical research of this system, it is necessary to set the certain boundary conditions. It 
can be made as follows. Two conditions, necessary for the beginning of inflationary stage, 
were specified above. First of these conditions can be written as the following restriction 



on value of inflaton and dilaton fields in the beginning of inflation [11]: 

e 2b *m 2 o 2 = 1. (4.5) 

The fulfilment of the second condition can be checked directly at numerical research of 
the set of equations ( |4.4| ). Thus evolution of the universe will pass in two stages: first, 
the stage of contraction will occurs from a point of the beginning of inflation up to a turn 
point. Then, after the bouncing off, the model will come on a stage of expansion. 

So, having use the condition (4.5) and, as well as in a case of closed model, setting 



some values of momenta q = qo and v = v$ in a point of the beginning of inflation we can 
carry out numerical research of the set of characteristics. Its results are presented in Figs. 




Figure 3: The dependence <j>{a) and a(a) for b = 0.2 (left panel) and for b = 1 (right panel). Here 
ai is a point of above the barrier reflection, a; - point of the beginning of inflation. 



5. Probabilities of tunneling and above the barrier reflection of wave func- 
tion 

In the considered above models, evolution of scalar fields was investigated during tunneling 
and above the barrier reflection. However, it is also interesting to find probabilities of these 
processes. It appears that it is possible to make in analytical form because of one feature 
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of the considered sets of characteristics. It is clear from the numerical analysis of Eqs. 
J4.2| ) in case of tunneling that action S is large at small values of momenta q and v. Thus 
permeability of a barrier will be exponentially small in accordance with the formula 



D = exp(-25). 



(5.1) 



Therefore the basic contribution to permeability will be brought by characteristics with 
large momenta for which the action is small. In this case it is possible to simplify essentially 
the set of characteristics (^^) because at large momenta's values their practically does not 
vary during tunneling j^, |17f| . Then Eqs. (4.2) will become: 

d<p q da p— bc Pi 

~ ~da~~~'a 2 F 

v = const, 



v 



da 



where q, v » 1, F 



a 2 F'' 
q = const, 
dS_ _W_ 
~da~ ~ ~F' 

A/ a and A 2 = q 2 + v 2 . Then it is easy to find 



(5.2) 



— In — . 

A a 



a 



v 
bq 



bg/A 



and 



A 



1 , a* 2^/4 
--In— + v 1 -(al 
4 ao 4 



m 2 v 2 2 1 6 



(5.3) 



(5.4) 



In these expressions ao is a point of an entrance of wave function under a barrier, a* is 
a point of an output. Last one is determined from the condition W(a*, (j)*, <r*) = and, 
finally, is a function of momenta q and v. 

As it has already been mentioned above, these formulas are true only at large values 
of q and v. However, in this case action S is small that does not allow to use the usual 



formula for permeability of a barrier (5J) which is applicable only at large S. Therefore 
in this case it is necessary to use the formula of generalized WKB [p!q| : 

-25) 



D 



exp( 



(5.5) 



l + exp(-25)' 

which is used at any S. Then, using Eqs. (|5.4|), ( |5.5| ) and averaging on all characteristics, 



(5.6) 



we can obtain the following expression for permeability of a barrier: 

e W (-2S(q,v)) 



D 



lim 

■ ■ q v 



dqdv 



o o 



l + exp(-2,S (q,v))' 



In a case of flat model the above the barrier reflection will occurs near to a maximum 
of superpotential ( 2.1l| ) where it is small on absolute value in comparison with the term 
(q 2 + v 2 )/a 2 from Eqs. ( f4.4|) . Then we shall have again the same expressions (5^) for 
fields. For calculation of reflectivity we use the known formula fl9||: 



R 



exp 



-41m 



(12 



day/-W(a) 



(5.7) 
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where a% is an arbitrary point on a real axis and 02 is a turn point calculated from the 
condition W = 0. In our case, using Eqs. ( |5.3| ) and ( p. 1 1 ) for k = 0, we have: 



W 



1 



2, ,2 



and 



4a 2 b 2 q 2 

1/3 



f bq Y /A 
° 2 = I J 6XP + 

Then the reflectivity, averaged on all momenta, will be: 



(5.8) 



(5.9) 







111 I 7T im V 

R = lim / / dqdv exp — tt -i / ~r - 



b q 



(5.10) 




The obtained formulas (5.6) and ( 5.10) ) 
allow us to find the probabilities at vari- 
ous values of parameter b (Fig. ||). As we 
see from figure, in case of the closed model 
permeability is the larger when b is larger 
(i.e. when the barrier is more thin), that is 
quite natural. On the other hand, as was 
noted above, than b is smaller, the infla- 
ton field on an output from under a barrier 
is larger and the conditions for inflation are 
better. So, certain compromise is necessary 
at a choice of parameter b for maintenance 
both enough long inflation and high proba- 
bility of tunneling. 

In a case of flat model the similar sit- 
uation occurs. Here, then b is larger, decreasing of the inflaton field a is faster on the 
expansion stage (see Figs. @). Then the conditions for inflation at identical values of the 
inflaton field in the point of above the barrier reflection a<i in a case of large b will be worse 
than at small b. At that, the probability of above the barrier reflection is increasing at 
large 6, as in a case of closed model. 

I am grateful to H.-J. Schmidt and H. Kleinert for discussions. This work was supported 
by ISTC project KR-677. 
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Figure 4: The dependence of permeability D 
and reflectivity R from b. 
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